Frequency dependent transport through a spin chain 
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Motivated by potential applications in spintronics, we study frequency dependent spin transport 
in nonitinerant one-dimensional spin chains. We propose a system that behaves as a capacitor for 
the spin degree of freedom. It consists of a spin chain with two impurities a distance d apart. We 
find that at low energy (frequency) the impurities flow to strong coupling, thereby effectively cutting 
the chain into three parts, with the middle island containing a discrete number of spin excitations. 
At finite frequency spin transport through the system increases. We find a strong dependence of the 
finite frequency characteristics both on the anisotropy of the spin chain and the applied magnetic 
field. We propose a method to measure the finite-frequency conductance in this system. 

PACS numbers: 75.10.Jm, 75.10.Pq, 75.30.Ds, 75.76.-fj 
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I. INTRODUCTION 

One of the main issues in modern electronics is that as 
devices get ever smaller the removal of waste energy gen- 
erated by Joule heating (which is in turn caused by the 
scattering of the information carrying conduction elec- 
trons) becomes problematic. A possible solution to this 
problem of excess heating is offered by spintronics^— in 
nonitinerant materials. Since charge transport is absent 
in such materials, the energy that is dissipated in the 
transport of a single information carrier (a spin) in a 
nonitinerant system^ is much lower than in both itiner- 
ant electronic- and spintronic devices^* 6 - Furthermore, it 
has been theoretically shown that the energy required for 
a single logic operation is low in a spin based device.— 

Clearly, energy considerations are not the only factor 
that play a role in the determination of the feasibility 
of spintronics in nonitinerant materials. Other impor- 
tant criteria are the operating speed of potential devices, 
and the possibility of creation, control, and read out of 
the signals used in these systems. Indeed, it has been 
proposed that spin-based devices in nonitinerant systems 
could operate at high frequencies.™ Furthermore, there 
have been various proposals to create (by means of spin 
pumping 8 and the spin Seebeck effect 9 ) and read out (by 
means of the inverse spin Hall effec t 10 ' 1 ) these pure spin 
currents in nonitinerant systems. Logic schemes that are 
based on interference of spin waves in insulating mag- 
nets have been proposed] 12 : 13 In these schemes, the dis- 
persion relation of the spin waves is controlled either by 
magnetic fields, magnetic textures,— electric fields due to 
the induced Aharonov-Casher phased or in multiferroic 
materials) 16 ' 17 or spin orbit coupling. 18 

In this work we focus on transport of spin excitations in 
antiferromagnetic nonitinerant one-dimensional quantum 
systems (spin chains). In such systems charge transport 
is absent, and the elementary excitations that transport 
magnetization are spin waves. Of the possible candidates 
for such systems that exist we specifically mention here 
CS2C0CI4, which has a large anisotropy in the exchange 
interaction,— and SrCu02, in which transport was re- 
cently shown to be ballistic over distances on the order 



of 1 ^m (see Ref. |2fl). 

Inspired by possible applications in spintronics, we set 
out to design the spin analogues in insulating magnets of 
the different components that are used in modern elec- 
tronics. Recently, we have shown that systems containing 
spin chains can be used as a spin resistance^ and a spin 
diode. 21 In this work, we propose a device based on spin 
chains that mimics the working of a capacitor, but for 
the spin degree of freedom. That is, a device that has 
a spin conductance that is zero for applied DC driving 
fields and increases with frequency for applied AC driv- 
ing fields. Traditional capacitors have a wide range of 
applications in electronics, most notably in frequency se- 
lection, noise reduction, and temporary energy storage. 

The device that we propose consists of a spin chain 
with anisotropy A y - in the exchange interaction in the z 
direction. We introduce two impurities, that each model 
the replacement of a single atom, in the spin chain. A 
renormalization group argument shows that the impuri- 
ties flow to strong coupling at low energies, or, equiva- 
lently, low frequencies. At low frequencies the chain is 
therefor in an insulating state, with an island containing 
a discrete number of spin excitations in the middle. We 
show that at higher frequency the spin conductance be- 
comes nonzero; the system undergoes a transition to a 
conducting state. We show that the point at which this 
transition occurs can be tuned by an external magnetic 
field and the bulk anisotropy in the exchange interac- 
tion. Transport through such double barrier systems has 
received attention recently* 2 ^ - — but these studies were 
restricted to the itinerant electronic case. Here we focus 
on the nonitinerant magnetic case with an emphasis on 
applications. 

In Section [TT] we introduce our system and the tech- 
niques we use to analyze the system. Explicit expressions 
for the finite frequency spin conductance of the system 
under consideration are derived in Section Mil This fi- 
nite frequency spin conductance contains all the required 
information to show that our system behaves as a spin 
capacitor. In Section IIVI we present numerical results, 
conclusions are given in Section fVl 
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FIG. 1. (a) Schematic view of the nonitinerant spin sys- 
tem. The 1-D spin chain of finite length L is adiabatically 
connected (meaning that the length of the transition region 
L t 3> 4a) to two spin reservoirs, (b) Illustration of the work- 
ing principle behind the spin capacitor. Strong impurities pin 
the field <j>(f), thereby blocking the system for spin transport. 



II. SYSTEM AND MODEL 

A spin chain of length L (see Fig. [T]), adiabatically 
connected to two three-dimensional spin reservoirs, can 
be described by the anisotropic XXZ-Heisenberg Hamil- 
tonian 



H = Y, Jij [SfSf + S?S] + + ////,,• ^ BiSf. 

(«> « 

w 

The summation is over nearest neighbors, and we avoid 
double counting of the bonds. Here Sf is the a- 
component of the spin operator at position Fj. Jy de- 
notes the exchange interactions between the two nearest 
neighbors at f\ and fj. We assume an, in general spa- 
tially varying, anisotropy Ay > in the exchange cou- 
pling. This anisotropy typically originates either from 
a Dzyaloshinksii-Moriya (DM) or dipole-dipole interac- 
tion between the different spins.— The last term in the 
Hamiltonian describes the Zeeman coupling between the 
external time-dependent magnetic field, Bi(t) = Bi(t)e z , 
and the spin operator. We assume the magnetic field to 
be of the form B t (t) = B - [1 - 9{-L/2)] ABcosut. 

For our spin capacitor, we need to introduce two im- 
purities in the spin chain. We propose to use an impurity 
of a specific type, a single instance of which, located at 



is described by its anisotropy 



A, 



A5ij-i (S iiio + 8i,i +i) = A.j 



(2) 



Such an impurity describes the replacement of a single 
ion in the spin chain by another ion with a different 
anisotropy. We choose this form because it contains all 
the important physics for a general impurity. With two 
impurities, located at = id/2, and allowing for a con- 
stant anisotropy A' (which we assume to be smaller than 



1), the anisotropy becomes 



Ay — A_ d / 2 



A 



d/2 



A'. 



(3) 



It is well knowr* 2 ^ that in the one-dimensional spin 
chain we can map the Heisenberg Hamiltonian on a 
fermionic Hamiltonian by means of the Jordan- Wigner 
transformation. Assuming that g/isB -^i J and Ay <C 1 
we can take the continuum limit in the resulting fermionic 
Hamiltonian, which we can then linearize around the 
Fermi wave vector. In the process we need to intro- 
duce left- and right going fermions V'l/^tr)- We can 
then perform a bosonization procedur e) 29 ' 30 using the 
density field 4>(f) and its conjugate field 9(f), which sat- 
isfy [cj)(x),d X '9(x')} = itt5(x — x 1 ). Using the fermionic 
creation operator for r = L/R- going particles 



4(x) = 



_ e ~ie r k F x e i[t r <j)(x)-9(x)] 



(4) 



where a is the distance between two neighboring spins, 
and e r = ±1, we arrive at the following real-time action 
describing the system 



S[0] = S , oM + S , ,M + 5bM. 



(•5) 



Here Sq [4>] is the quadratic action which originates from 
the in-plane terms and the forward scattering terms in 
the out-of-plane part of the Heisenberg Hamiltonian. It 
is given by 



2ttK 



d 2 r 



u{d x 4>{f)f--{d t 4>{r)f 
u 



(G) 



Here u is the propagation velocity of the bosonic excita- 
tions and K is the Luttinger Liquid interaction parameter 
which can be determined from the parameters J and A' 
of the spin chain from the Bethe-Ansatzi^i We denote 
f = (x,t). Sb [0] describes the Zeeman coupling between 
magnetic field and the z component of the spin operator 



9VB 



A 2 rcj){r)d x B{r). 



(7) 



Lastly, S z [4>] contains the most relevant contributions 
coming from the two impurities. We have shown previ- 
ously that, for an impurity of the form Eq. @ centered 
at x , the most relevant terms are given by2i 

Sx [<t>] = ^2 J dt [Xl cos20(xo,i) + A^ o sm2^(x ,t)l8) 

Here, the couplings A^o^ are of the order J A. Further- 
more, a — Kg{iBBo/{tujJ c ) is the doping away from half- 
filling due to the magnetic field Bq. uj c is the UV-cutoff 
of the theory, which is on the order of J/ h. We have thus 
S z [(f>] = £-(2/2 M + Sd/2 [</>]• Note that we must assume 
here that 2kpd ^ tt, since otherwise the two impurities 
would cancel each other. To lowest order the backscatter- 
ing terms flow under renormalization as 1 — K . Since we 
work in the regime 1/2 < K < 1 (which corresponds to 
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< A' < 1) the impurities flow to strong coupling for low 
energies. In this regime, the spin chain is effectively cut 
into three parts. The middle island (—d/2 < x < d/2) 
contains a discrete number of excitations, and tunneling 
through the island is only possible by co-tunneling pro- 
cesses, in which two excitations hop on and off the middle 
island simultaneously. These processes are strongly sup- 
pressed at low temperatures. For finite frequency, pro- 
cesses similar to photon-assisted tunneling become possi- 
ble, so that we expect the spin current through the spin 
chain to increase as a function of frequency. 

From Eq. ([8]) it is seen that the strength of the impuri- 
ties is proportional to the filling a. Hence we can control 
the finite frequency characteristics of the system by ap- 
plying a constant magnetic field Bq to our system. Most 
strikingly, this field allows us to control the frequency 
at which the system will become conducting in the spin 
sector. Since the Fermi wave vector also depends on Bq, 
this field also controls the number of spins that are on the 
island for a given width d of the island. The number of 
spins relative to half-filling can be estimated as (d/a)a. 

Lastly, we can model the excitations in the two reser- 
voirs by the quadratic Luttinger Liquid action as well. 
We will assume two reservoirs that are described by the 
the Heisenberg Hamiltonian with A = 0, so that the ex- 
citations in the reservoirs are free, and we have K r = 1. 



III. FINITE-FREQUENCY SPIN 
CONDUCTANCE 

As long as we consider only magnetic fields that 
point along the quantization axis, the spin current I s (f) 
through the system can be determined from the conser- 
vation of spin, and is given by 



d t (j>{r). 



(9) 



The nonlocal differential finite frequency conductiv- 
ity a(x,x',tu) is defined as the Fourier transform 
<j(x,x',uj) — J dte lu)t a(x, x' , t) of the real-time linear re- 
sponse of the spin current to an infinitesimal change in 
the time-dependent magnetic field gradient 



a(x, x' , t — t') — 



s(W)) 

6[d*B(i*)Y 



(10) 



In circuit theory, the quantity which is of interest is the 
spin conductance G(w), which we define for our system 
as 



G{u>) =Re [a(-L/2,L/2,u)\. 



(11) 



With this convention a positive magnetic field difference 
AB yields a positive spin current for the free system. 
To calculate the transport properties we use the Keldysh 
formalism.— We assume that at t = — oo the system is 
described by the quadratic action Sq[4>] + Sb [4>] and that 



the perturbation S z [<fi] is turned on adiabatically. We cal- 
culate the conductivity to lowest non-zero order in S z [4>] ■ 
In this description the conductivity is given by 



u(x, x ,t — t) = i 



g^B S ( SZ[J] 
7T 6[d x ,B(r>)] \6J(f) 



(12) 



j=o 



Here, Z[J] is the partition function of the system [see Eq. 
(|Alj) ] . and J(r) a generating functional. We give a short 
explanation of the evaluation of Eq. (|T2|) in Appendix [A] 
and summarize the important results here. We can show 
that the conductivity consists of two contributions 

a(x, x',u) = <r {x, x',u>) + <Jbs(x, x',u). (13) 

The first contribution to the conductivity describes the 
system in the absence of impurities. It gives the exact 
spin current for the free system, since in the absence of 
impurities the current is linear in magnetic field gradient. 
The free conductivity is given by 



cr (x, x', w ) = - 



( 9HB 
\iry/K 



(iu))iG$(x,x',u>). (14) 



The retarded Green's function for x, x' in the chain or 
at the boundaries can be obtained from analytic contin- 
uation of the Matsubara Green's function for the inho- 
mogeneous Luttinger Liquid models If we assume in- 
finitesimal dissipation in the leads, the retarded Green's 
function is given by 

n Ri ' \ ttK ( „ I 
G Q (x, x , (jj) 



2uj \ 



sy— 1 g — iCb 



E 



=±i 



D irLu(x-\-x') 



— irLd\x~x' | _|_ , 



(15) 



Here 7 = (1 — /s)/(l + k) is the reflection coefficient 
for Andreev-like (in the sense that for some processes 
the transmission coefficient can be larger than 1) den- 
sity wave scattering at the chain- reservoir interface, and 
k = K/K r . The dimensionless frequency Co — uj/ujl, 
where ujl = u/L is the energy scale determined by the 
finite length of the spin chain. Also x = x/L. In the limit 
u) — > we have a(x, x', lu) = K r (g^B) 2 /h as expected. 

The quantity aBs(x,x' ,uj) contains corrections to the 
conductivity coming from the impurities. For the situ- 
ation with two impurities, it depends on the amplitude 
and both the absolute and relative position of the impu- 
rities, and is given by 



a B s(x 7 x',uj) = 



(—) E (16) 



As it turns out, the terms with £ = £' reduce to the 
known result for a single impurity , 34 i 35 the terms with 
£ 7^ £' come from interference between the impurities 
located at —d/2 and d/2. Here a^^i(x,x' ,w) is given by 

<r (x,£,w) 



a^(x,x',u}) = -- 



■[a (e,x',u)iK XZ,?,u>) 



-(T {Z,x',u)iK{Z,S',0j\. 



(17) 
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K(£,£',w) denotes the retarded backscattering current- 
backscattering current correlator, which is defined by 

/>oo 

iK({,g,u) = dre^ ([I BS ^,r),I BS (e,0)})- (18) 
Jo 

The backscattering current operator is defined as 

gfj, B 5S Xo [<f>] 



ttH 6<j)(xo,t) 



[Xi cos20(x o ,i) - X a xo sin2</)(xo,a]9.) 



In taking the functional derivative, (j>(xo, t) is regarded as 
a function of t only. Some lengthy but straightforward 
algebra shows that we can rewrite Eq. (IT8l) as 



iK{^ uj) = 2 (^) 2 [A|A|/ + A|A|,] dre^ 



,2C (f,£',T) _ „2C Q (£',£ ,-t) 



(20) 



where Co(x, x', r) is the regularized ^-correlator at zero 
temperature for the free Luttinger Liquid action on a 
system of finite length L, given by 



C (x,x',t) = 
K 



2 1 

n £ even 



Win 



n 2 + a 2 



7i G odd 



l lln 



(n + x — J') 2 + a 2 (ir + l) 2 

[(n + 2i) 2 + a 2 ] [(n + 2£') 2 + a 2 
[(n + x + x') 2 + a 2 (if + l) 2 ] 2 



21) 



Here f = uj c t, a = ujl/uj c , and w c is the UV-cutoff of the 
theory, which is on the order of J/h. Finally, we have to 
replace the bare couplings A in Eq. (|20|) by the renor- 
malized couplings, given by (max[wi,w] /w c ) _1+A A. In 
the next Section we will evaluate Eq. (fT6|) numerically. 



IV. RESULTS 

We have shown in Sec. IIIII that the spin conductance 
G(uj) = Re [er(— L/2, L/2, uj)] of our system consists of 
two parts, the free conductance Gq(uj) and the backscat- 
tering conductance Gbs( w )- These are of opposite sign 
for low frequencies. Fig. [2^ shows a plot of log 10 |Gbs ( w ) I 
(Gbs(w) is here in units of (g/is) 2 /h) as a function of 
log 10 (ai/cj c ) for different values of the Luttinger liquid 
parameter K. Different values of K correspond to dif- 
ferent values of the constant anisotropy A' in the spin 
chain (see Eq. Q). It is seen that Gbs(w) goes to 
zero at high frequencies. In the high-frequency regime 
the spin conductance is therefore determined by the free 
conductance Go(w) only, and hence the system is con- 
ducting. For lower frequencies the backscattering con- 
ductance Gbs(w) grows, so that the total spin conduc- 
tance decreases. Since our calculation of Gbs(w) is per- 
turbative in the strength of the impurities, we cannot 
trust our calculations at frequencies uj < uj* , where uj* is 
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FIG. 2. (top) Plot of log 10 |GbsH| ( in units of (jCu) 2 /h) 
as a function of log 10 (o;/a; c ) for different values of the Lut- 
tinger liquid interaction parameter K (increasing from top to 
bottom). Used parameters are B = 70 mT, J = 10 2 K, 
L = 10 6 a, d = 10 4 a, K r = 1, and A = 0.3. The black dots 
denote the frequencies uj* where |Go(w*)| = |Gbs(^*)|- Our 
perturbative results cannot be used for u < uj* . (bottom) 
Idem for different values of Bo (decreasing from top to bot- 
tom). Used parameters are K = 0.7, J = 10 2 K, L = 10 6 a, 
d — 10 4 a, K r = 1, and A = 0.3. The oscillatory behavior is 
explained in the caption of Fig. [3] and the associated part in 
the text. 



the frequency such that |Gbs(w*)| = |Go(w*)|. Indeed, 
instead of the behavior shown in Fig. [5] the total spin 
conductance G{uj) is approximately zero at frequencies 
below uj* . Therefor, uj* is approximately the frequency 
at which the system undergoes the transition from insu- 
lating to conducting behavior. 

Interestingly, we see that the backscattering conduc- 
tance is strongly dependent on the Luttinger liquid pa- 
rameter K. If we focus again on the conductance at the 
frequency uj* , we see that an increase from K = 0.6 to 
K = 0.7 leads to a decrease of the frequency at which 
the system switches between conducting and insulating 
behavior of approximately 2 orders of magnitude. We 
note that we have uj c « J/h, which for J = 10 2 K is ap- 
proximately uj c s» 2.1 • 10 12 s _1 . The previously discussed 
increase in uj* is therefore approximately from 7.4 * 10° 
rad s _1 to 6.7 * 10 7 rad s _1 . 

We see from Fig. [5] that the frequency uj* strongly de- 
pends on the applied magnetic field Bq as well. This is 
caused by the fact that the strength of the impurity is 
proportional to Bq, see Eq. ([5]). This strong dependence 
is convenient for possible applications, since it would al- 
low for an externally tunable capacitor. For the parame- 
ters in Fig. [2] we see that by increasing Bq from Bq = 70 
mT to Bq = 700 mT, we can increase uj* from 7.4 * 10 s 
rad s _1 to 9.5 * 10 7 rad s _1 . 

In Fig. [3] we show the dependence of Gbs(w) on the 
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FIG. 3. Backscattering conductance Gbs(^) as a function 
of ui/ujc for different values of the impurity separation d (in 
units of a). Used values are B = 1.4 T, J = 10 2 K, L = 10 4 a, 
A' = 0.7, K r — 1, and A = 0.2. It is seen that, independent of 
d, the conductance oscillates with period 2-kujl (we estimate 
oil ~ a/L here, the values for the free system). The values of 
d only influences the shape of the oscillations, not the period. 



pling at low energies. By including two such impurities in 
the spin chain we got a device which is insulating at zero 
frequency driving field AB, and has a spin conductance 
that grows in magnitude under an increase in frequency 
of the driving field. We have studied the influence of the 
anisotropy in the exchange interaction A and the exter- 
nal magnetic field Bo, and have found that both have 
a strong influence on the finite frequency characteristics 
of the system. We have proposed a way to measure the 
effects. 
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distance d between the two impurities. Regardless of d, 
it is seen that the conductance oscillates with frequency 
2hlol (see caption). These oscillations are caused by in- 
terference of the bosonic excitations due to Andreev-like 
reflections at the boundaries between reservoir and spin 
chain. It has been shown that, in the presence of a dou- 
ble barrier structure in the chain, additional oscillations 
with frequency determined by the parameters d/L and 
(L — d)/ (2L) are visible in the current at finite driving 
field. 26 These oscillations do not show up in the conduc- 
tance at zero AB; however, they are expected show up in 
the conductance at finite AB. Here, d changes the shape 
of the oscillations, but not the frequency. All the calcu- 
lations here were done at T — 0. Generally, the effect of 
finite temperature is to wash out the fluctuations. 

The finite frequency spin conductance could be mea- 
sured by measuring the spin accumulation in a system 
consisting of a spin diode of the form proposed in Ref. 
l2ll in series with the spin capacitor proposed here, be- 
tween two spin reservoirs. We can apply the driving field 
AB(t) to the left reservoir, and measure the spin accu- 
mulation in the right reservoir. In this setup, the spin 
accumulation is zero when the spin capacitor is in the 
insulating state. As was shown in Ref. 21, if we con- 
sider 10 4 parallel spin chains with K — 0.6 in a magnetic 
field Bq = 750 mT, and amplitude of the driving field 
AB = 43 mT, we have a spin accumulation of approx- 
imately 10 12 magnons per seconds if the capacitor is in 
the conducting state. 



Appendix A: Partition function 

In this Section we discuss the partition function of the 
system described in Sec. |TTJ For simplicity, we consider 
a system that is described by the action So[<p] + Sb[4>] + 
At cos [20(0, t)\. The partition function of the 
system described in Sec. [TT| can be inferred from the 
results for this action in a straightforward manner. The 
partition function can be written in vectorized form as 



Z[J{r)} 



_D0exp 



d 2 rdV x 



exp 

where we have defined 
${r) 



cj>{r) ■ G (r, ?)<!>(?) - 2iS(r- r , )(j)(r) ■ QJ^r 1 ) 



dt V cos [2^(0, t)} 



(Al) 



4>+{r) 
cj>-{r) 

J W = ^( J(r1/2 



G (r>S) 



J_ 1-1 
y/2 V 1 1 

G +(r,r<) G^-(r^) 



(A2) 



V. CONCLUSIONS 

In this work we have proposed a spin capacitor in a sys- 
tem consisting of one-dimensional nonitinerant spin chain 
adiabatically connected to two spin reservoirs. The spin 
chain is required to have an anisotropy in the exchange 
interaction. We have shown that the replacement of a 
single atom in the spin chain leads to an local backscat- 
tering term in the Hamiltonian that flows to strong cou- 



Here Gq' v (r, f ) is the contour-ordered Green's function 



So 



(A3) 



Here, 77, 77' = ± denotes whether a field is located on 
the positive or negative Keldysh contour. Eq. (TT2"j) can 
be evaluated by first making the linear transformation 
4>'{r) — 4>{r) —if Ar l Go(r, r')Q T J(r l ) and then perform- 
ing the functional derivatives^ 4 - 
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